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Abstract 

The vector space V'' of the eigenfunctions of the Laplacian on the 
three sphere 5*^, corresponding to the same eigenvalue Xk — —k {k + 
2), has dimension (fc + 1)^. After recalling the standard bases for V*, 
we introduce a new basis B3, constructed from the reductions to 5*^ of 
peculiar homogeneous harmonic polynomia involving null vectors. We give 
the transformation laws between this basis and the usual hyper-spherical 
harmonics. 

Thanks to the quaternionic representations of and SO(4), we are 
able to write explicitely the transformation properties of B3, and thus 
of any eigenmode, under an arbitrary rotation of S0(4). This offers the 
possibility to select those functions of V* which remain invariant under 
a chosen rotation of S0(4). When the rotation is an holonomy transfor- 
mation of a spherical space S^/F, this gives a method to calculates the 
eigenmodes of S'^/F, which remains an open probleme in general. We 
illustrate our method by (re-)deriving the eigenmodes of lens and prism 
space. In a forthcoming paper, we present the derivation for dodecahedral 
space. 

1 Introduction 

The eigenvalues of the Laplacian A of 5^ are of the form Xk = —k (fc + 2), 
where k G IN"^. For a given value of k, they span the eigenspace V*^ of 
dimension {k+ 1)^. This vector space constitutes the {k+ 1)^ dimensional 
irreductible representation of SO (4), the isometry group of S^. 

There are two commonly used bases (hereafter Bl and B2) for V'' 
which generalize in some sense (see below) the usual spherical harmonics 
Yem for the two-sphere. The functions of these bases have a friendly be- 
havior under some of the rotations of SO (4); this generalizes the property 
of the Yim to be eigenfunctions of the angular momentum operator in IR"^. 
However, these functions show no peculiar properties under the general 
rotation of S0(4). 

Excepted for some cases (lens and prism spaces, see below), the search 
for the eigenmodes of the spherical spaces of the form /F remains an 
open problem. Since those are eigenmodes of which remain invariant 
under the rotations of F, it is clear that this search requires an under- 
standing of the rotation properties of the basis functions under SO (4). 

The task of this paper is to examine the rotation properties of the 
eigenfunctions of V*^ , as a preparatory work for the search for eigenfunc- 
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tions of /r (in particular for dodecahedral space). This will be done 
through the introduction of a new basis B3 of V*' (in the case k even), 
for which the rotation properties can be explicitely calculated: following 
a new procedure (that was already applied to in ^) we generate a 
system of coherent states on V*. We extract from it a basis B3 of V*, 
which seems to have been ignored in the literature, and presents original 
properties. Each function "l>jj of this basis B3 is defined as [the reduction 
to of] an homogeneous harmonic polynomial in IR'*, which takes the 
very simple form {X ■ N)'' . Here, the dot product extends the Euclidean 
[scalar] dot product of IR'' to its complexification C*, and A'^ is a null vector 
of C**, that we specify below. After defining these functions, we show that 
they form a basis of V'' , and we give the explicit transformation formulae 
between B2 and B3. 

The properties of the basis B3 differ from those of the two other bases, 
and make it more convenient for particular applications. In particular, it 
is possible to calculate explicitely its rotation properties, under an arbi- 
trary rotation of S0(4), by using their quaternionic representation (section 
0. This allows to find those functions which remain invariant under an 
arbitrary rotation. In section |1J we apply these result to (re-)derive the 
eigenmodes of lens and prism space. 

2 Harmonic functions 

A function / on is an eigenmode [of the Laplacian] if it verifies A/ = 
A/. It is known that eigenvalues are of the form At — —k {k + 2), k £ IN"*". 
The corresponding eigenfunctions generate the eigen [vector] space V'°, of 
dimension (fc+ 1)^, which realizes an irreducible unitary representation of 
the group SO (4). 
First basis 

I call Bl the most widely used basis for V'' provided by the hyper- 
spherical harmonics 

Bl = {Vu^ oc Yim), I = l..fc, m = -LI. (1) 

It generalizes the usual spherical harmonics Yim on the sphere. In fact, 
it can be shown (pp, |5] p. 240, [2]) that a basis of this type exists on any 
sphere 5". Moreover, [5| show that the Bl basis for 5" is " naturally 
generated " by the Bl basis for 5""^. In this sense, the Bl basis for 5*^ 
is generated by the usual spherical harmonics Yim. on the 2-sphere S'^ . 

The generation process involves harmonic polynomials constructed 
from null complex vectors (see below). The basis Bl is in fact based on 
the reduction of the representation of S0(4) to representations of S0(3): 
each yhim is an eigenfunction of an S0(3) subgroup of SO(4) which leaves 
a selected point of 5*^ invariant. This make these functions useful when 
one considers the action of that peculiar SO(3) subgroup. But they show 
no simple behaviour under a general rotation. We will no more use this 
basis. 

Second basis 

By group theoretical arguments, IT" construct a different ON basis of 
1/*, which is specific to S'^: 

B2= (Tfe;„,,„J, 7711, 7712 = -fc/2...fc/2, (2) 

where mi and 7712 vary independently by entire increments (and, thus, 
take entire or semi-entire values according to the parity of fc). In the 
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spirit of the construction refered above, B2 may be seen as generated 
from a different choice of spherical harmonics on S^. The bases Bl and B2 
appear respectively adapted to the systems of hyperspherical and toroidal 
(see below) coordinates to describe S"^. 

The formula (27) of m, reduced to the three-sphere, shows that the 
elements of this basis take a very convenient form if we use toroidal coordi- 
nates (as they are called by |7|) on the three sphere S^: {x, 8, 4>) spanning 
the range < X ^ ""/S, < 8 < 2n and < < 27r. They are conve- 
niently defined (see [3 for a more complete description) from an isometric 
embedding of 5*^ in IR* (as the hypersurface x £ IR''; | x \— 1): 

a;" = r cos x cos 6 

x^ = r sin x cos 4> 

x^ = r sin x sin (f> 

x^ = r cos X sin 9 

where (a;^), /i — 0,1,2,3, is a point of IR''. As shown in jj), these 
coordinates appear naturally associated to some isometries. 

Very simple manipulations show that, with these coordinates, each 
eigenfunction of B2 takes the form: 

^fc;mi ,m2 (-^) — ^fc;mi,m2(x) ^ ^ ' (^) 

where the tfc;mi,m2(x) ^re polynomials in cosx and sinx and we wrote, 
for simplification, £ = mi -I- m2, m = m2 — mi. 

To have a convenient expression, we report this formula in the har- 
monic equation expressed in coordinates X) 8, (j). This leads to a second 
order differential equation (cf. equ 15 of [?])• The solution is proportional 
to a Jacobi polynomial: tk;mi,m2ix) cos^ x sin"" x ^d"'^(cos2x), d = 
k/2 — m2. Thus, we have the final expression for the basis B2 

Tk:„^i,m2{X) = Ck-^mum2 [cos X e^f [siu X e"*]™ P^'"'^'[cos(2x)], (4) 

with Cfe;™i,„2 = ^^^^ \/ [t^/llZfy. tk/l-'Zly. from normalization re- 
quirements (the variation ranges of mi anf m2 imply that the quantities 
under factorial sign are entire and positive). Note also the useful propor- 
tionality relations: 

COS X sm X P k-i-m) (cos 2x) oc cos x sm X P l-i+rL) (cos 2x) 

2 2 

OC COS X sm X P k+\-.r. (cos 2x). 

2 

The term C C = e"* e™"^ in ^ defines the rotation properties of 
Tfc;mi,m2 under a specific subgroup of S0(4). This properties generalizes 
the properties of the spherical harmonics on the two-sphere , to be 
eigenfunctions of the rotation operator Px . This advantage has been used 
by |7] to calculate (from a slightly different basis) the eigenmodes of lens or 
prism spaces (see below sectionllj. However, the Tk-^mi,m^ have no simple 
rotation properties under the general rotation of S0(4). This motivates 
the search for a different basis of V''. 

Note that the basis functions Tk-^mi,m2 have also been introduced in 
(p. 253), with their expression in Jacobi Polynomials. Note also that 
they are the complex counterparts of those proposed by 17] (their equ. 19) 
to find the eigenmodes of lens and prism spaces. The variation range of 
the indices mi,m2 here (equ. |5J is equivalent to their condition 

\l\ + \jn\<k, l + m = k, mod (2), (5) 

through the correspondence I = mi + m2, m = m2 — mi. 
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2.1 Complex null vectors 

A complex vector Z = {Z^ , , Z'^ , Z^) is an element of C*. We extend 
the Euclidean scalar product in IR'* to the complex (non Hermitian) inner 
product Z ■ Z' = J2^ {Z'Y, fi = 0,l, 2, 3. A null vector is defined 
as having zero norm N ■ N = N^N'^ = (in which case, it may be 
considered as a point on the isotropic cone in C*). It is well known that 
polynomials of the form {X ■ A)*^, homogeneous of degree k, are harmonic 
if and only if A is a null vector. This results from 

Ao(A • A)*^ = ^ {X-Nf i Y^iNf. A^) j (A • A)*"^ = 0, 

where Ao is the Laplacian of IR*. Thus, the restrictions of such polyno- 
mials belong to V'°. As we mentioned above, peculiar null vectors have 
been used in [5| and to generate the bases Bl and B2. 

To construct a third basis B3, let us first define a family of null vectors 

N{a,b) = {cosa, i sinfe, i cost, sin a), (6) 

indexed by two angles a and b describing the unit circle (they define 
coherent states in IR**). 

The polynomial [X ■ N{a,b)]'' is harmonic and, thus, can be decom- 
posed onto the basis B2. It is easy to check that, like the scalar product 
X ■ N{a,b), this polynomial depends on a and b only through the com- 
binations e''*""' and e^^'^'^''\ with their conjugates. This implies that its 
decomposition on B2 takes the form 



ib(m2 —m\ ) 



[A- A(a,b)]'= = ft;™i,™2 rfc;mi,m.(A) e-''"(™i+'"^) e 

(7) 

where the coefficients Pk;mi,m2 do not depend on a,b. Now we intend to 
find a basis of V'' in the form of such polynomials. 



2.2 An new basis 
2.2.1 Roots of unity 

To do so, we consider the (fc -I- l)**^ complex roots of unity which are the 
powers of 

p = e = cos Q -f I sm a, a = . (8) 

fc -I- 1 

We recall the fundamental property, which will be widely used thereafter: 

k 

Y P"' = (k + l) <5f (9) 

n = 

where the equallity in the Dirac must be taken mod fc -|- 1. 
In a given frame, we consider the family of null vectors 

Nij = N{Ia, Ja) = {cos la, i sin Jq, i cosja, sin 7a), I,J — 0..k 

(10) 

and we define the functions $jj : $Jj(A) = (A • Nu)^ . We report such 
a function in equ.Q to obtain its development in the basis B2. Then we 
muhiply both terms by p^(™i+™2)-./(m2-mi) Making the summations 
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over /, J (each varying from to fc), and using @, we obtain, in the case 
where k is even (that we assume hereafter): 

1 

T, — I(mi+m2)-J(m2-mi) 

where ^7^-^^ — -ffc;mi,m2 ^fc;mi,m2* 

This gives the decomposition of any rfc;mi,m2 (and thus, of any har- 
monic function) as a sum of the (fc + 1)'^ polynomials providing the 
new basis of V*^ : 

B3 = {^''u), I,J = 0..k (fceven). (12) 

The coefficients Pk-mi.m2 involved in the transformation are calculated 
in Appendix A. We obtain easily the reciprocal formula expressing the 
change of basis: 

fc/2 

$?J= Yl 7'fe;mi,m2 + (13) 

mi ,m2 — — fc/2 

3 Rotations in IR^ 

3.1 Matrix representations 

The isometries of 5*^ are the rotations in the embedding space JR*. In the 
usual matrix representation, a rotation is represented by a 4*4 orthogonal 
matrix g G S0(4), acting on the 4-vector (x^) by matrix product. 

In the complex matrix representation, a point (vector) of IR^ is repre- 
sented by the 2*2 complex matrix 

; W = x° + ix^, Z = x^+ix^£ C. 

A rotation g is represented by two complex 2*2 matrices {Gl,Gr), so 
that its action takes the form X > Gl X Gr (matrix product). The two 
matrices Gl and Gr belong to SU(2). Since SU(2) identifies with S^, any 
matrix Gl or Gr is of the same form than the matrix X above. Since 
SU(2) is also the set of unit norm quaternions, there is a quaternionic 
representation for the action of SO(4). 

3.2 Quaternionic notations 

Let us note j^, ^ — 0, 1, 2, 3 the basis of quaternions (the correspond to 
the usual l,i,j,k but we do not use this notation here). We have jo = 1- 
A general quaternion is q — q^j^i — g'^ -I- q^ji (with summation convention; 
the index i takes the values 1,2,3; the index takes the values 0,1,2,3). 
Its quaternionic conjugate is q = — q^ji. The scalar product is q\ ■ q2 = 
((?i q2 + q2 gi)/2, giving the quaternionic norm \ q \'^= ^ = I]p('7'^)^- 

We represent any point x = (a;^) of IR* by the quaternion q^ = x'^ ji^. 
The points of the (unit) sphere correspond to units quaternions, | g p = 
1. Hereafter, all quaternions will be unitary (if no otherwise indicated). It 
is easy to see that, using the coordinates above, a point of 5*^ is represented 
by the quaternion cos X C + sin x C Ji > where we define dotted quantities, 
like ^ = cos 9 + js sin 0, = cos (j) + ji sin 4>, as the quaternionic analogs 



X ; 



W_ iZ 
iZ w 
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of the complex numbers = cos^ + i sinO and ^ = cosiji + i sm0, i.e., 
with the imaginary i replaced by the quaternion js. 

In quaternionic notation, the rotation g : x i-^ gx is represented by a 
pair of unit quaternions {Ql,Qr) such that q^: qgx = Ql qx Qr- 

Complex quaternions, null quaternions 

The null vectors A'' introduced above do not belong to IR** but to 
C*. Thus, they cannot be represented by quaternions, but by complex 
quaternions. Those are defined exactly like the usual quaternions, but 
with complex rather than real coefficients. Note that the pure imaginary i 
does not coincide with any of the j^, but commutes with all of them. Also, 
complex conjugation (star) and quaternionic conjugation (bar) must be 
carefully distinguished. Then it is easy to see that the (null) vectors Njj 
defined above correspond to the complex quaternion nij = + i j2 p'^ ■ 
Note that | mj 0. 

In quaternionic notations, the basis functions are expressed as 

^ij(x) = (Nij ■ X) =< njj ■ qx > = ^ ^ ) ' 

Quaternionic notations will help us to check how the basis functions are 
transformed by the rotations of SO (4). 

3.3 Rotations of functions 

To any rotation g, is associated its action Rj on functions: Rg : / i-^ 
Rg/; Rg/(a::) = f{gx). Let us apply this action to the basis functions; 

TZg$ij{x) = ^ij{gx) nij ■ (Ql q^ Qr) . (15) 

We consider a function on 5*^ also as a functions on the set of unit quater- 
nions {qx is the unit quaternion associated to the point x of 5*^). On the 
other hand, we may develop this function on the basis: 

fc 

^,9ij=Y. GYAg) (16) 

The coefficients G^j j{g) of the development, that we intend to calculate, 
completely encode the action of the rotation g on the basis B3, and thus 
on V''. 

To proceed , we introduce three auxiliary complex quaternions: 

Q = 1 + i J3, /? = ji - i j2 = (1 - i js) Ji and 5 = -ji - i j2. 

They have zero norm and obey the properties < a ■ nij >— p^ , 
< a ■ nij >= p^' , < f3 ■ nij >= p'' , < 5 ■ njj >= Let us now 

estimate the relation H16^ for the specific quaternion a + Ra + SP + T5, 
with R,S,T arbitrary real numbers: 

(A+RA'+SB+TV)'' = J2g'iA9) <{p'+Rp~'+S p'+T p-^)\ (17) 

where A =< Ql ct Qr- nu >, A' =< Ql a. Qr ■ mj >, B =< Ql (3 Qr ■ 
nij >, V =< Ql S Qr-uij > characterize the rotation. (Note that these 
quantities depend on / and J). 

We develop and identify the powers of the exponents R, S, T: 

A" A'"-" d'^-p-'- = J2 G7j{g) p'^^"-"'' pj{2'-fe+p). 



6 



This holds for < q < p, < r < k — p, < p < k. After definition of 
the new indices A = q + r, B = q — r + k — p, which both vary from to 
k, the previous equation takes the form 

A'B) \BV ) \ A 



AB 
A'V 



i(A+B-k)+i(A-B) 



This holds for any value of A,B,p. A consequence is that AA' = BV, 
which can be checked directly. Finally, 

ij 

withW=(f),V=(#). 

Taking into account the properties of the roots of unity, this equation 
has the solution 



{A') 



k 

(fc + l)2 ^ ^ 

^ ' A,B=0 



i(A+S-fc) -j{A-B) yl 



(18) 



When a rotation is specified, there is no difficulty to estimate the associ- 
ated values of A! , U, V, and thus of these coefficients which completely 
encode the transformation properties of the basis functions of under 
S0(4). 

In the next section, we apply these results to rederive the eigenmodes 
of Lens or Prism space. In the next paper |^, we take for g the generators 
of r, the group of holonomies of the dodecahedral space. This will allow 
the selection of the invariant functions, which constitute its eigenmodes. 

4 Lens and Prism space 

The eigen modes for Lens and Prism space have been found by 0. Here 
we derive them again for illustration of our method. 



4.1 Lens space 



An holonomy transformation of a lens space takes the form, in complex 
notation. 



Gl = 







, Gr = 







Its action on a vector of IR takes the form 



X 



W iZ 
iZ w 



Gl X Gr — 



G. 
(19) 

(20) 



In this simple case, W = + ix^, Z = x^ + ix^ are transformed into 
W e^^^ and Z e^^^ respectively. This corresponds to the quaternionic 
notation 



}l = wi W2, Qr = wi/w2, Wi = cos(i/'i/2) + sin(7/),/2). 



(21) 



The rotation is expressed in the simplest way in the toroidal coordi- 
nates, since it acts as S i— > ^ -I- t/ji , (f) (j) + il)2- From the expression 
of the basis functions (B2) , it result their transformation law : 
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This leads directly to the invariance condition Itpi + mtp2 ~ mod 2tt. 
Using the standard notation for a lens space L{p, q), namely 



= 2-n/p, i>2 = 2n q/p, 
we are led to the conclusion: 



the eigenmodes of lens space L{p, q) are all linear combinations of Tk-mi,m2, 
where the underlining means that the indices verify the condition 
mi + m2 + q{ m2 — mi ) — 0, modulo(p). 



4.2 Prism space 

The two generators are single action rotations {Gr = 0). The first gener- 
ator, analog to the lens case above, with i^i = tp2 ~ 2n/2P, provides the 
first condition £ + m = 0, mod 2P which takes the form 

rn2 = mod P. (22) 

This implies that k must be even. 

r — i 

The second generator has the complex matrix form G = Gl ~ . g 

which corresponds to the quaternion Ql ~ Q = — ji- Easy calculations 
lead to ^ = —p'', A' — p~'\ B — , V — —p^^. Reporting in 11811 gives 

G% = 7^7^i: (-+^-^+^^) (-1)^. (23) 

y'^^'-' A,B=0 

This formula, together with those expressing the change of basis between 
B2 and B3, allow to return to the rotation properties of the basis B2 which 
take the simple form: 

R. n- (_'i\m.2+k/2 q- /"O/l"! 

It results immediately that the G-invariant functions are combinations of 

-'fe;mi,m2 V / -'fc;m;^, — m2- 

Finally, 



the eigenfunctions of the Prism space are combinations of 
Tfc;„i,„2 + (-l)IIi2+fc/2 t;,^^ Vmi; k even. 



According to the parity of fc/2, the functions 7fc;mi,o are included or not, 

from which simple counting give the multiplicity as 

{k + 1) (1 + [fc/2P]), for k even ([...] means entire value), 

(fc + 1) [fc/2P], for k odd, in accordance with [J. 



5 Conclusion 

We have shown that V*^, the space of eigenfunctions of the Laplacian of 
with a given eigenvalue \k (k even) admits a new basis B3. In contrary 
to standard bases (Bl and B2) which show specific rotation properties 
under selected subgroups of SO (4), it is possible to calculate explicitely 
the rotation properties of B3 under any rotation of the group S0(4), 
as well as to calculate the functions invariant under this rotation. This 
opens the door to the calculation of eigenmodes of spherical space. The 
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eigenfunctions of lens and prism spaces had been calculated by |7|, by 
using a basis related to B2 (its real, rather than complex, version). We 
rederived them to illustrate the properties of the bases. 

In a subsequent paper we apply these results to the search of 
the eigenfunctions of the dodecahedral space S^/F, where F — Dp is 
the binary dihedral group of order 4P. Those functions, still presently 
unknown, are the eigenfunctions of which remain invariant under the 
elements of F. 

5.1 Appendix A 

Let us evaluate the function 

k 

ZL{X) = P""'" *'^(^) (25) 
/./=o 

1 1 1 

cosx (Cp^'^ + +sinx (?P'' - ^) , 

where we defined = e'* and ^ = e**^. After development of the power 
with the binomial coefficients, the sum becomes 

E P"'"^' E ( ' ) (CP-^ + [sinx (Cp-^ - 7^)]^ (26) 

IJ p=0 ^ ^ / 

Let us write the identities 

(^p-^ + 1 )'=-P = ( ^"7^) p-^(2'-+''-'=-^), (27) 

^P r=0 ^ ^ 



= 2 2^ p 



^P q=0 V / 

that we insert into I26I I. After summing over I, J, and rearranging the 
terms, we obtain: 

7 (X. _ o-fc u V (-1)^"" (cosx)''-^^+'" (sinx)^^-" 

This formula results from the fact that, through JSJ, the summations over 
I, J imply p = 2q — m and 2r — £ + k + m — 2q, that we have reported. 
The range of the summation over q is defined by the conditions 

< £ + k + m - 2q < 2k + 2m - Aq < 2k, 0<q<2q-m<k. (30) 

Rearrangements of the previous formula, inserting u = cos(2x) ~ 
2 cos^ X ~ 1 = 1 ~ 2 sin^ Xi ^^Sid to 

(m + d)! (£ + d)! 

E ("^^') (rft-) (1+^)^-1)^- 
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where we have defined i = ^+"^-^ — g and d = ^^'^^"^ . Verification shows 
that the range defined as above gives exactly the development formula for 
the Jacobi polynomial. The comparison with gives the coefficient 

P, ^ rUil 

^^^2 - mi)! (fc/2 + mi)! (fc + 1)^ Ck-m,,r^, 
_ 2-'' TT k\ (fc+1)-^/^ 

^ V(fc/2 + m2)\ (fe/2 - m2)!(fc/2 + rm)! (fc/2 - rm)! 
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